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Dispersive estimates for rational symbols and local well-posedness 
of the nonzero energy NV equation 

Anna Kazeykina* and Claudio Munoz^ 


Abstract 

We consider the Cauchy problem for the two-dimensional Novikov-Veselov equation in- 
tegrable via the inverse scattering problem for the Schrodinger operator with fixed negative 
energy. The associated linear equation is characterized by a rational symbol which is not a 
polynomial, except when the energy parameter is zero. With the help of a complex analysis 
point of view of the problem, we establish uniform decay estimates for the linear solution with 
gain of almost one derivative, and we use this result together with Fourier decomposition 
methods and X s,b spaces to prove local well-posedness in H s , s > 1. 


1 Introduction. Main results 

In this paper we consider the initial value problem associated to the Novikov-Veselov (NV) 
equation, in two dimensions, with a fixed energy parameter E € M (28j : 

dtv = 8(d z + dg)v + 2 d z (vw) + 2 d z (vw) — 2 E(d z w + d z w), (1.1) 

d z w = —3d z v. (1.2) 

Here v = v(t,x,y), w = w(t,x,y), and 

z = x + iy = (x, y) € C. 

The complex partial derivative operators d z and d z are defined as usual: 

d z := i( d x - id y ), d z := i( d x + id y ), (A x>y = 4d z d z = 4:d z d z ). 

Formally, one has w = — 3d^~ 1 d z v, where d^ 1 d z can be defined via the Fourier transform V, 
acting on L 2 (M 2 ), £ = (£i,£ 2 ) € M 2 \{(0,0)}: 

m"D 2 /](£i,6) = (|^|)^[/](b,6)- (1-3) 

In addition, note that the right hand side in (11.11) has the structure of ‘a real part’. Consequently, 
one looks for v real-valued: 

v(t, x,y)sM, w(t, x, y) € C, 
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so that the first equation can be written as follows: 

dtv = 4Re{d 2 (4d|u — 3vd^ 1 d z v + 3 Ed^ 1 d z v)}. 

For the sake of completeness, we write equations (O-dLlD in terms of real-valued coordinates 
(x,y). If we identify w = w\ + 11 V 2 with the vector held w = ( W\,W 2 ), with w\,w 2 real-valued, 
then m becomes 

dyW\ + d x u >2 = 3 d y v, dyU >2 — d x w\ = 3 d x v, 

and m reads 

dtv = 2 \d x (d 2 v — 3 d^v) + V.(vw) — EV.w] . (1.4) 

See Appendix[A]for a proof of this fact. However, for most of this paper we will work with the z-z 
formulation (ll.ll) - (ll.2D . which presents several advantages when computing precise estimates. 

The NV equations are important for several reasons. First of all, (ED-D is a completely 
integrable model in two dimensions, in the spirit of the Kandomtsev-Petviashvili (KP) and the 
Davey-Stewartson (DS) equations [19] , However, unlike the two last models, the origin of NV is 
not physical, but mathematical [28]: it arises as the integrable equation obtained by assuming 
that the associated scattering problem corresponds to the standard, stationary Schrodinger 
equation in two dimensions with fixed energy E : 

(-& x ,y + v{t,x,y) - E)ip = 0. (1.5) 

In that sense, NV is the closest model that generalizes (in a mathematical form) the famous 
Korteweg-de Vries (KdV) equation to the two dimensional case, and it has interesting connec¬ 
tions with the scattering problem for bounded or localized potentials, and the theory of inverse 
problems for sufficiently smooth potentials, see e.g. [8] and references therein for more details. 
Moreover, if v depends only on x , we have z = z = x, d z = d z = ^d x and 

w = — 3v, dt.v = 2d x (d x v — 3v 2 + 3 Ev), 

which is a KdV equation^ 

Another interesting connection between NV and other integrable models is the following. It 
can be formally shown (see jGr2j for example and Appendix [B] for more details) that when the 
parameter of energy E tends to ±00 the NV equations, rescaled appropriately, become in the 
limit KP equations (KPI corresponding to E —>• +00 and KPII corresponding to E —» — 00 ). 
It was also shown formally in [Gr2j that there is a a corresponding convergence of scattering 
problems for the KP and NV equations. In order to obtain a rigorous proof of such asymptotic 
regime, a suitable well-posedness theory for NV equations is necessary, with explicit bounds of 
the solution and its lifespan depending on E. 

The Novikov-Veselov equation was first obtained in an implicit form by S.V. Manakov in 
[22] , It has the following operator representation 

d t L=[L,A] + BL, (1.6) 


where 

L := -A X}V + v{t, x, y) - E, A := -8(d)? + df) - 2 (wd z + wd z ), 


and 

B := 2 (d z w + d z w), 

1 Note that in this case u(t, x ) := —v(—t/ 2, x — 3Et) solves the standard KdV equation ut + u XX x + 6uu x = 0. 
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with w given by (11.21) and [•, •] denoting the commutator. Equation (1 1 .6|) is the compatibility 
condition of 

Lp = 0, tp t = Ap. 

Representation (11.61) is called Manakov triple representation for (11.11) and can be considered 
as a generalization of the Lax pair representation for KdV (see [20]) to the (2 + l)-dimensional 
case. 

From the integrability of (1 1.1 D - (1 1.2 1) it follows that NV formally has an infinite number of 
conserved quantities. Here we write the first three that have a correct meaning (each integral is 
taken on C = K 2 ): the L * 1 integral 

J v{t)dxdy = J v(0)dxdy, (real-valued), (1.7) 

the “mass” 

M[v\{t) := J vw{t)dxdy 

= J vw(0)dxdy, (complex-valued, two identities), (1.8) 


and the “energy”: 




J 6 d z wd z v + Eu? 


vw 


(' t)dxdy 


i7[u](0), (complex-valued, two identities). 


(1.9) 


We see that iL 1 (M 2 ) is the energy space, for if v £ H 1 , then d z v and d~v are in L 2 , and w £ H 1 . 
It is not difficult to check that we also hav(0 


v,w € L 3 (M 2 ). 


Compared with other dispersive models coming from Physical theories, NV equations lack of 
signed conserved quantities which control the long time dynamics. For this reason, a good 
understanding of some particular explicit solutions is essential. These regimes will strongly 
depend on the sign of the energy parameter E : we will see that E < 0 corresponds to a sort 
of defocusing case, whilst E > 0 is the focusing regime. The case E = 0 is different in several 
aspects. 

Indeed, the literature [19] describes three different regimes associated to (fLTD-(fL2D, depend¬ 
ing on the value of the energy E: the case where E = 0 (usually denoted as NVq), and the cases 
where E > 0 (NV+) or E < 0 (. NV- ). 

In the case E = 0, using direct and inverse scattering techniques, Perry [29] has been able to 
show existence of a solution for a certain set of initial data. Angelopoulous [1] showed that the 
Cauchy problem for (11. 1 1) - (1 1.2 D is locally well-posed for v(t = 0) £ H s ( R 2 ), s > by following 
the ideas of Molinet and Pilod [24] for the Zakharov-Kuznetsov equation. Below L 2 (M 2 ) he 
showed [lj that the flow has some ill-posed behavior. One may ask if by using (11.81) or (11.91) the 
well-posedness result may be extended to a global one. The answer is surprisingly no (another 

2 Indeed, this fact is a consequence of the Gagliardo-Nirenberg inequality 

IMIlp <C[|Vn||2 2 H^“, ot £ [0,1), p = ———• (1.10) 

1 — a 
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fact that reveals that NV does not enjoy standard physical properties). For any a, c, d € M such 
that a + c(x 3 + y 3 ) + d(x 2 + y 2 ) 2 > 0 everywhere, the function 

v(t, x , y) = -2A XjV log (a - 24ct + c(x 3 + y 3 ) + d(x 2 + y 2 ) 2 ) 

solves NVq, decays like r^ 3 at infinity (r = \f x 2 + y 2 ), and it blows up at finite time (see also 
[33]). Note additionally that if v(t,x,y) is solution of NV, then for all A > 0, 

v\(t, x, y) := X 2 v(X 3 t, Ax, Ay), and w\(t, x, y) := X 2 w(X 3 t, Ax, Ay) (1.11) 


satisfy (I1.1D - (I1.2I) with 


E x := EX 2 


( 1 . 12 ) 


In the case E = 0 the equation does not change and we have another solution of the same 
equation]! Using this scaling symmetry, we can construct a solution that blows up at every time 
of the form t\ , A ^ 0, at the same points as before. Moreover, since Ed 1 or even L 2 are subcritical 
regularities, one realizes that there are arbitrarily small solutions that blow up in finite time, 
unlike in standard dispersive models. In addition, the time of existence will not depend on the 
size of the solution only. 


In the case of negative energy, NV- is in some sense reminiscent of the KPII equation studied 
by Bourgain in [5]. It has been shown in m via the inverse scattering techniques that in this 
case the equation does not possess sufficiently regular soliton solutions decaying as r -3 ~ £ , e > 0 
(compare with results of [3], [4], [2] on the absence of sufficiently localized solutions of KPII). 
The inverse scattering theory also allowed to obtain the large time asymptotics (uniform in 
space) for sufficiently regular solutions of NV- with nonsingular scattering data (in particular, 
sufficiently small solutions have nonsingular scattering data) (see [13] and compare with an 
analogous uniform large-time asymptotic result for KPII in [18)1. In this paper, we prove the 
following result. 


Theorem 1.1. NV- is locally well-posed in H S (M'?) for any s > The lifespan of solution is 
proportional to the size of the energy \E\. 


The precise version of this Theorem is formulated as Thereom 15.11 below. Here we present 
briefly the ideas of the proof of this result. 

In the proof of Theorem ll.ll we follow the ideas of Bourgain, extended and simplified by many 
authors; in particular we use the formalism employed by Molinet and Pilod in [24] fot the study of 
low-regularity solution to the Zakharov-Kuznetsov (ZK) equation in several geometries in 2 and 
3 dimensions. Note that in [24| the authors deal with symbols whose resonance functions seem 
a priori intractable via Bourgain’s ideas (their method was developed in [1] for NV in the case 
E = 0). It is essential in their work that Molinet and Pilod use a sharp L i Strichartz estimate 
with gain of derivatives obtained by Carbery, Kenig and Ziesler [7] for polynomial symbols like 
the one in ZK. In our case, we do not have a polynomial symbol, but a fractional one, much in 
the spirit of KP equations. We overcome the difficulty of showing a Carbery-Kenig-Ziesler type 
estimate by proving almost sharp smoothing estimates in two dimensions for the linear part of 
the equation, in the spirit of Kenig, Ponce and Vega [16], and Saut [30J, a task that requires 
some fine Fourier analysis coming from the improvement of previous results showed by the first 
author mm in the framework of the Inverse Scattering Transformation approach. It turns 
out that we gain almost one derivative in L°° with decay slighty slower than 1 /t, which suffices 
for us to close standard Bourgain’s bilinear estimates for s > 

3 Therefore, the critical space for E = 0 is H~ 1 . 
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In addition to these estimates, we perform a standard Fourier localization between low 
and high frequencies. We recall that the resonance function that appears when dealing with 
the interaction of low-high to high frequencies is treated by estimating the zero level set via 
reasonable lower bounds on the partial derivatives of the resonance function. Such estimates are 
simple to establish, but they carry a loss of accuracy (probably 1/2 of derivative) that could be 
avoided by dealing directly with the resonance function as Bourgain did in [5]; however, such a 
task substantially more difficult given the complexity of the linear NV symbol (see (14. ip for more 
details). However, it is worth to mention that the nonzero energy NV symbol has some useful 
boundedness properties near the origin, unlike standard KP equations. Finally, the high-high 
to high frequencies interactions are treated using the almost sharp smoothing estimate for the 
linear symbol, which leads to a Strichartz estimate with 1/4 gain of derivative in L 4 . 

We finish this introductory section with some words about the global existence problem. 
Although the equation is in nature well-behaved for high regularity initial data, the problem of 
global existence is far from trivial because of the lack of an evident sign in the real and imaginary 
part of the conservation laws, see (11.81) and m- Even worse, some nice initial data may have 
bad behavior in terms of conserved quantities: for example any smooth, rapidly decaying radial 
data has radial Fourier transform and consequently it mass m is identically zero: 

/ /•OO /*27T 

vw= / e~ 2l9 \v(r)\ 2 r drdd = 0. 

Jo Jo 

This implies that in principle the mass cannot be used to bound the L 2 norm of the solution, as 
usually performed in other physical models. A similar conclusion holds for the two first terms 
in the Hamiltonian (11.91) . 

It turns out that in the case of NV equations, such a problem is deeply related with the 
behavior of scattering solutions of the associated Schrodinger operator (11.51) . The results of 
Novikov m imply global existence for initial data with suitable spectral properties, via the 
inverse scattering transform. Very recently, Perry [29j showed that in the case of zero energy, if 
the initial data are of conductivity type, in addition to other decay assumptions, then IVVo has 
a global weak solution that is classical if the initial data are in the Schwartz class. Recall that 
a potential vq € Cq°(M 2 ) is of conductivity type if the equation 

(-A x ,y + vo)tp = 0 

admits a unique positive solution such that </? = 1 in a neighborhood of infinity. This definition 
comes from the theory of Inverse Problems, see [6] GHEE] for more details on this interesting 
connection. 

Finally, some words about the case E > 0. The methods of this paper are not entirely suitable 
for the case E > 0. The main obstacle is the fact that the smoothing Strichartz estimate is 
not longer sufficient to control the large-to-large frequencies regime, and an approach of the 
type Carbery-Kenig-Ziesler is probably needed. We hope to consider the remaining case in a 
forthcoming publication. 

Acknowledgments. The second author would like to thank the Laboratoire de Mathe- 
matiques d’Orsay for their kind hospitality during past years, and Didier Pilod for explaining his 
paper [211 an d giving several useful comments about a preliminary manuscript. The authors also 
thank Jean-Claude Saut for the initial motivation of this study and for many useful discussions. 

Notation. In the following text notation “A < R” means that there exists a constant c 
(depending only on a) such that A < cB. 
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2 Smoothing estimates for negative energies 

The aim of this section is to estimate the integral 

I(t, u; E) := J |£| Q+i/ y^(^ ;E )dRe£cflm£, (2.1) 

c 


where 


and 


Note that 



S(u,£',E) — (£ 3 + 1 3 ) ^1 — + u 0- 

p(0 = (£ 3 + f) (i-jpf) 


is the symbol of the linear part of the NV equation. 


( 2 . 2 ) 


Lemma 2.1 (Dispersion estimate for negative energy). Let 0 < a < 1 and (3 € R. Then for all 
t > 0 and for any fixed e(a) > 0 small, the following estimate is valid: 


I I(t,u-,-l)\ < 


(1 + 1/31) 

^(o+3)/4—e ’ 


uniformly on u € C. Recall that the implicit constant depends 


on a only. 


(2.3) 


Remark 2.1. Note that we obtain better estimates than in the case of the Zakharov-Kuznetsov’s 
oscillatory integral, see e.g. [21]. Compared to the results in [21] . when a ~ 1 we gain almost 
one derivative instead of half a derivative. 

Remark 2.2. By making a suitable change of variables, we easily prove a dispersion estimate in 
the case of arbitrary negative energy E < 0. We have 


\nt,u-E)\ < 


(l + I^DI-El ^- 1 )/ 8 - 36 / 2 

t(«+3)/4—£ 


(2.4) 


with constants independent of E. Indeed, given (12.11) with phase (12.21) . we are reduced to the 
case E = — 1 by performing the change of variables £ = l-E] 1 / 2 ^, to = | E\ 3 ^ 2 t, and = |-E| 1 / 2 ^, 
which converts m into 


\E\ a / 2 \C 0 \ a e ito§ ^°’-^\E\dRe^ dlm£ 0 , u 0 := g. 


Using (12.31) we are lead to 


\I{t,u-E)\ 


(l + lfl)!^ 


t, 


q+3 

4 

0 


3(^+3) _3 

<1+ 10)1*1*? |B| * 1 


q+3 

t~ - £ 


< 


(1 + 


i -.-.I a — l 3 

I E\ 8 2 £ 

q+3 7 ; 

t~i — £ 


as desired. Note that in the formal limit E —>• 0, estimate (1231) becomes singular. 
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Some useful consequences of Lemma 12.11 are stated below. 


Corollary 2.2 (Smoothing and Strichartz estimates). Let U(t) = U(t ; E) be the associated NV- 
linear group, namely for £ = fi + if, 2 , 

U{t)v 0 := f e «S(«>^V 0 
Jr 2 

(c/. (12.11) and (12.21) 1. T/ien /or any 0 < a < 1, f3 € [0,1], we have the smoothing decay estimate 


|||0,r^(i)« O || I £ <a,/3 


and the Strichartz estimate with gain of almost one-half derivative: 


\\\d z \ a ^ 2 U(t)v 0 \\m L P < a>/3 


L l,y' 


(2.5) 


( 2 . 6 ) 


with p = -SLp, i + = 1, and | = /3(^ -e). 

Proof. The proofs are by now standard, but we prove it for the sake of completeness. See e.g. 
ehhb] and references therein for detailed proofs. 

First, we prove m- This is just a consequence of the standard complex interpolation 
theorem. We have from (12.41) and Young’s inequality, 


\m a+ip u(t)v 0 \\ L - < 


tH^- £ 


ToIIli, 


for any 0 < a < 1 and f3 € R. We interpolate against the trivial estimate 


\\\d z fU(t)v 0 \\ L2 = \\v 0 \\ L ,, 


I r <? tP ? 


to get (|2.5D . 

Now we prove ». By duality, we are lead to prove that 

J f>{t,x,y)\d z \ als/ 2 U(t)vodtdxdy < \Ef^ i £ ) | 
for all <f € Cg°(M 3 ) and with | ^ = 1. However, 

J (f(t, x,y)\d z \ a ^ 2 U{t)vQdtdxdy = Jv o J \d z \ a/ 3 ^ 2 U(—t)(j)(t,x,y)dt dxdy 

\d z \ a P ! 2 u {—t)4>{t)dt 


< 


L 2 


L 2 


Now, 


|a z | a ^ /2 Z7(-t)</(i)dt 


L 2 


< 


[ [ [ \d z \ a, 3 , 2 U(-t')(j)(t',x,y)\d z \ a P / 2 U(-t)(t)(t,x,y)dxdydt'dt. 

[ [ [ \d z \ a ^U{t' — t)cj>(t',x,y)(f(t, x, y)dt'dtdxdy 

J x,yJtJ t' 

[ | d z \ afi u{t'--W)dt' 

Jt' 


\m\ tc'jv' ■ 

t q t p m v-.y 

n t 
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On the other hand, we have from ([2.5ft and the Hardy-Littlewood-Sobolev’s inequality, 


)4>{t')dt' 

< 

f 


T q T p 

Jt' 


[ \d z \^u{t'--W) 

Jt r 


T p 

- L/ x,y 


dt' 


Li 


< 

r^j 


< 

r^j 


|£|( 3 (V-~! e ) 


L 


IWOII 


tP 

x ’ y dt' 


t' I t-tf {*¥-*) 


T q ' T p ' 

- LJ + J - / rr.. ?/ 


Another set of standard but useful estimates is the following. 


Corollary 2.3. We have 

ll C/ ( i )-ollL U <l^|- 1 / 16 + l|t/(^o|l x o^ + , 

and 

\\\dz\ 1/4 ~U(t)v 0 \\ L 4 < \E\-° + \\U(t)v 0 \\ i- 

The constant in the last inequality becomes singular as \ approaches 
Proof. Put a = 0, (5 = (| — e) 1 = | + in (12.61) . we get 


II U(t)v 0 \\i 


< 


t,x,y 


\m\v 0 \\ 


L l,y' 


with r = ^ 


i-fs 
!-| £ 


T + and 7 = -jg 


From here we have (see e.g. 


1 — 12 s \ 

1 -fe J 

Lemma 3.3]) 


j_+ 

14 


□ 


(2.7) 

( 2 . 8 ) 


II^IU u , < 1 ^ 11 ^ 11 ^,^, 

By interpolation with the trivial estimate for ||u||r2 , using the interpolation parameter 0 = 

t,x,y 

|(1 — |) = | and pg 1 = (1 — 9)pQ 1 + Qpf 1 , 0 < po < Pi, we obtain (12.71) . 

Now we deal with (12.81) . Taking a = 1 — 16e and f3 = ^77 = in (12.61) we have p = q = 
4 ^_|~ = 4 + . We also have tjo/3 = = \ and 

m^ir^oii^^iEi-^iiuoii^, 


with r = 4 


izif 

1—5s ’ 


from where 


ll|a,l r t/(f)»ollL U , < |Er 0+ ||(7(i)ooll x ,,... 

We interpolate with Lf x y with 6 = 1 — |e to obtain 

II|S,I i/ 4 ^(‘)»oIIl* <|£|- 0+ ||[/(tKH 

c ) a: )y _a ’ 2 


□ 


In the next Section we prove Lemma 12.11 
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3 Proof of Lemma 12.1 


3.1 Restatement of the problem 

This section is devoted to the proof of the pointwise dispersion estimate for nonzero energies. 
Compared to previous works (see e.g. mm), we estimate the whole two dimensional integral. 
We start by performing the following change of variable^: 


£ = -i(A- = ), f = z ( A - 


(3.1) 


Denote 


Note that 


¥>(A) : = -*( A - = ) =£. 


^:{A€C : |A| > 1} —> {£ € C 2 : £ ^ 0} 


is a bijective, smooth map. Note also that the Jacobian determinant of this transformation 
satisfies 

£>(£,£) 1 |A| 4 — 1 , , 

' = 1 - jyjl = Mtu— (3-2) 


Note also that 


D( A, A) 


\t\ a = 


a - i 

A 


|A | 4 

|AA-1|° 

|A| q 


(3.3) 


Finally we remark the almost trivial fact (but very important in the following computations): 
one has 

A = _f 
A £’ 

From (12.21) . (13.11) and (13.41) . we have 


(3.4) 


iS{u,£, — 1 ) — * ( (£ 3 + £ 3 ) ( 1 + TTjy ) + + £ u ) 


_3 

l ^ 2 


— * ( £ 3 T + 3£ = + 3£ - ) + - ( (A — = I w — (A —- I u 






A 1 


— i ( i\ A — 3——h 3^7 — ) — i\ A — 3——I- 3-^2 — ^3 


A ' A 2 A 3 
+i ( 2>i ( A — =| = — 3i ( A — 


A 1 


1\ A 


+- 


X) A 

A-= )u- (A— 


X A 


A 


u 


4 Recall that we work with E = —1. In the case of general E < 0, the right change of variables is given by 




A 


Note that this formulation breaks down when E approaches the zero energy level. 
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Rearranging similar terms, we obtain 



This last expression will be important for us. We denote it 

S(». A) - (a* + - A 3 - £) + i ((a - i) 8 - (a - i) «) ■ (3.5) 


Compared with 5(n, £, — 1), in the new expression 5(«, A) t/ie terms in A and in A are decoupled. 
The advantage of this fact is that in terms of A the stationary points of the oscillatory integral 
I(t,u ; —1) satisfy an algebraic equation of only one complex variable (see (13.71) b while in terms 
of £ the stationary points are defined by a system of two coupled algebraic equations (of two 
real variables). 

Moreover, note that S(u, A) is always a pure imaginary number, since it is the difference 
between a complex number and its complex conjugate. 

Thus, replacing m , m and (1331) in (T 2 TTT) . we are brought to estimate the following 
integral 

J(t, <*;-!)= JJ ~ ~ — e‘ s( "- A >dReA<iImA, (3.6) 

C\Si(0) 

where Ri(0) is the open ball centered at zero, of radius one. We also remark that estimates of 
integrals of type (13.61) have been obtained in mm , see also previous works [ISIU5] concerning 
the KdV equation in one dimension. However, this time we need an improved estimate because 
we need a suitable gain of derivatives in order to close the iteration argument. 


3.2 Study of stationary points 

We will start by introducing the following parametrically defined sets of the complex plane (see 
Fig. [[]). Let 

U := {u € C : u = 6{2e~ iip + e 2iip ), ip € [0, 2 t r)}, 

and 

U := {u € C : u = 6t(2e"^ + e 2 ^), t € [0,1] tp € [0, 2 vr)}, 

be the (closed) region enclosed by the curve U. These sets will be essential to understand the 
stationary points of the phase function S(u, A). 

With these definitions in mind, let us describe the properties of the stationary points of the 
function S(u, A) defined in (13.51) . These points satisfy the equatiorU (here S\ stands for the 
partial derivative with respect to A) 


^ u u . 9 3 i 

5a= 2 "2W “ 3A + F = 0 ' 


5 Here the symbol = means that equality to zero is satisfied for stationary points only. 


(3.7) 
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Figure 1: The curve U and its enclosed region U in the complex plane. Note that the U (and its 
interior and exterior) is invariant with respect to the transformations and z H > ze 2lk7r / 3 , 


A: = 0,1,2. 


Additionally, the degenerate stationary points obey the equation 

Saa = J - 6 A - ^ = 0. (3.8) 

We denote ( = A 2 , and 

^ u u . 3 

Q {u,Q = -- — -X + -^- 

(Compare with (13.71) . 1 Clearly, for each £ = Q(u ), root of the function Q(u,C), there are two 
corresponding stationary points of S(u, A), given by A = ±y/C- Since Q(u,Q has only three 
roots counting multiplicity (say Aq(u), Af (u) and A|(u)), in terms of the variable A 2 , the function 
S\(u, A) can be represented in the following compact form 

S\{u, A) = — (A 2 - Ag(u))(A 2 - A?(u))(A 2 - A|(«)). (3.9) 

Concerning the behavior of the roots A j(u), defined in (|3.9I) . we have the following result, see 
El for a proofJl 

Lemma 3.1 (Description of stationary points, |lll). Assume that u £ C is a fixed parameter. 
Then the following are satisfied. 

1. If u = 18e — —, k = 0,1,2 (see the vertices ofU in Fig. \3.iA) . then 

A o(u) = Ai (u) = A 2 (u) = e~^, 

and S(u,X ) has two degenerate stationary points, corresponding to a third-order root of 

. . _ 2irik 

the function Q(u , (), (1 = e 3 • 

®Note that in m the parametrization of the set U was slightly different, here we give a more precise one. 
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2. Ifu&U ( i.e. u = 6 ( 2 e llf + e 2llf )) and u / 18e 3 ; for k = 0 , 1 , 2 , f/ien 


A 0 (u) = Ai («) = e^ 2 , A 2 (n) = . 

Thus S(u,X ) das two degenerate stationary points, corresponding to a second-order root of 
the function Q{u,Cf), Ci = e*^, and two non-degenerate stationary points corresponding to 
a first-order root, C 2 = e~ 2z,p . 

3. If u E int U, then 

Xi(u) = e llfi , and \{u) / A j(u) for i / j. 

In this case the stationary points of S(u,X) are non-degenerate and correspond to the 
square roots of the roots of the function Q(u,f) with absolute value equals 1. 

4 ■ Finally, if u G C\U, then 

A 0 (u) = (1 + u)e iip/2 , Ai (u) = e~ itp , X 2 {u) = (1 + a;)" 1 ^ 2 , 

for certain ip € R and u > 0. 

In this case the stationary points of the function S(u , A) are non-degenerate, and correspond 
to the roots of the function Q(u,() that can be expressed as Co = (1 + T )z l<f , Ci = , 

C 2 = (1 + r )~ 1 e i<p , and (1 + r) = (1 + w) 2 . 

Note that Lemma O describes the behavior of stationary points with respect to the curve 
U. Note that when estimating (13.61) . one needs to take into account that I contains a multiplier 
that vanishes on the unit circle and grows at infinity. The neighborhoods of |A| = 1 and infinity 
being special regions, we will need to regroup the four cases above in three different ones, as 
expressed below. 


3.3 Idea of the proof 

Consider the integral I = I(t,u, —1) in (13.61) expressed as follows: 


I = I(t, u, — 1) = J f(X)exp(tS(u,X))dReXdlmX, (3.10) 

C\Bi(0) 

where / is a given function that vanishes on dB i( 0 ), the boundary of the unit ball centered at 
the origin. 

In what follows we suppress the differential dReAdlmA unless it is explicitly specified. We 
also anticipate that / is given by the expression 


/(A) 


||A|2- 1 |°+tf(|A|4- 1 ) 

|A|“+ 4 +*/3 


(3.11) 


In order to estimate uniformly on u € C the large-tinre behavior of the integral /, in the following 
subsections we will use the following general scheme. 


1. Fix e > 0. Consider D e as the union of disks with radius e and centers at the stationary 
points of S(u, A) in the closed set C\Ri(0). 
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2. Represent I = I(t, u, —1) in (13.61) as the sum of integrals over D e \Bi(0) and C\(Ri(0)UL> e ): 


I = hnt + I ext, where 


lint — 


/(A) exp(f5(w, A)), and 


AA-BRO) 


(3.12) 


I ext — 


/(A) exp(tS(w, A)). 


€\(B 1 (0)UD e ) 


3. Find an estimate of the form 


| lint | < e P R\ 


(uniformly on u, t ), for some /3,7 > 0, and where R is the maximum of absolute values of 
the stationary points of S. 


4. Integrate I ex t by parts using the Stokes formula [] to obtain 

lext = I\ + R + H, 


with 


II ' 21. 


dD e \B 1 (0) 


h :=- 


/(A)exp(t5(u, A)) - 
S\(u, A) 

fx(X)exp(tS(u,X)) 


C\(£i(0)UDe) 


Sx(u,X) 


and 


h ■= 


C\(Bi(0)UD e ) 

5. For each R hnd an estimate of the form 

e s R K 

\h\ ~ —» 


/(A) exp(fS(u, A))5 aa(m, A) 
A)) 2 


(3.13) 

(3.14) 

(3.15) 

(3.16) 


<5, k > 0, 


using careful estimates on particular regions of the plane. 

6 . Choose e depending on R and t so that a) Rnt and I ex t are bounded with respect to R 
and b) I in t and I ex t decrease (as t —>• 00 ) with the same speed. 


3.4 Proof of the estimate for large t 

3 

We will suppose that t > e !-<* (see subsection 13.51 for the explanation of the choice of the 
constant). 

First, we introduce some notation. Note that for /(A) defined in (13.111) we have ( f\ = d\f) 

/a (A) = /a,i(A) + /a, 2 (A) + /a, 3 (A), 


where 

(a + i/3)||A| 2 — l| a+i/? - 1 A(|A| 4 - 1) 

/A,lv / • ma+4+tf 


(3.17) 
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Figure 2: Standard setting for Case 1, with one pair of roots on the circle |A| = 1, another 
pair inside the disk |A| < 1 and the last pair outside the disk |A| < 2. The set D e is the dark 
shadowed region. 


/A, 2(A) := 


2||A| 2 - 1|"+^A 


and 


/A, 3(A) := - 


|^|«+2+i/3 

{a + 4 + i/3) ||A| 2 - 1|“+^(|A| 4 - 1)A 


(3.18) 


(3.19) 


2|A| q + 6 +*/ 3 

In the following we consider three different cases, depending on the values of the parameter u: 

1. Case 1. Here u E C\U is such that |Ao| = 1 + w > 2 (see Lemma l3Tl item|4}. 

2. Case 2. Here bGU (corresponding to items [U [2] and [3] in Lemma 13.111 , and finally, 

3. Case 3. The region where u E C\U and |Aq| = 1 + u < 2 (see Lemma 13.11 item 13} . 


3.4.1 Case 1 


In this case the stationary points Ao and — Ao are separated from the unit circle, but ±Ai (u) lies 
on the unit circle. Fortunately, the third root A 2 (u) is uniformly outside the region of integration, 
since |Ao | = 1 + w > 2. 

We take D e to be the union of disks with radius e (to be chosen later) and with centers in 
points Ao, — Ao, Ai, —Ai (see Fig. ©. We also define 

R := |A 0 | = l+ w, (3.20) 

and we assume that e < R. 

We start by estimating Ii nt from (13.1211 . Note that for A E D e \B i(0) we easily have 

1 < |A| < |Aq| + e = R + e<R, 


and therefore the following estimates are valid: 


(iAi 2 -ir 

|A|" 


|A|“ < R° 


(3.21) 
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Figure 3: The set C\fi from Case 1 is the white open region. 


and 


|A| 4 —1 
|A| 4 


= 1 - 


|A| 4 - L 


(3.22) 


Thus, using the fact that S(u, A) in (13.51) is a pure imaginary quantity, and using (13.111) . (13.211) 
and (|3.22p . 


L int 


< 

/ 

l/(A)|| exp(tS(u, A))| 


JD e \B i(0) 


< 

f 

( A 2 — 1)"( A 4 — 1) 


Jd^BRO) 

|A|«+ 4 

< 

R a 1 

dlle A r/ImA. 


De\B 1(0) 


We conclude that 

\hnt\ < £ 2 R a - (3.23) 

We pass to the estimate of I ex t- We will start by splitting each of the integrals I 2 and I 3 into 
two pieces, the first one being an integral over the set of A at a distance greater than some fixed 
value from the stationary points, and the second one over the rest of the domain C\(.D £ U-Bi(0)). 
More precisely, we set 


C\n := {Ae C: |A| > 2 , |A — Aq| > 1 , |A + A 0 | > 1}, 


and consequently, 

n\(D e U Bi(0)) = {A € C : e < |A ± A 0 | < 1} 

u({A € C : 1 < |A| < 2} n{A £ C : |A± Ai| > e}). (3.24) 

(see Figures El El and 0]) and we denote /±, integrals with the same integrands as I 2 , I 3 
respectively and with the domain of integration being C\fi for the + sign and f l\(D e U i?i(0)) 


15 












for the — sign. More precisely, from (13.151) 


J+ := _ 


1 f /A(A)exp(tS(u, A)) 


C\f2 


In : =- 

2 t 


S\(u, A) 


/a(A) exp(tS(u, A)) 


f2\(D e USi(0)) 


5a(«,A) 


and from (13.161) . 


and 


/+ ■= I [ /( A ) exp(f5(n, A ))S X \(u, A) 


C\f2 


ir : = 


t 


n\(D e UBi(0)) 


(5a(«,A)) 2 


/(A) exp(t5(u, A))5 A a(«, A) 
(5 A (n,A)) 2 


(3.25) 


(3.26) 


(3.27) 


(3.28) 


We first treat integrals 7(7 (our goal is to obtain an estimate for these integrals indepen¬ 
dent of u). For that we use the following estimates valid on C\fl (see Fig. [3]): 


A 4 — 1 

< i 

A 2 — 1 

A 4 

rKJ ± ’> 

l A l 


|(A-Ai)(A + Ai)(A-A 2 )(A + A 2 )| 
|A| 4 


|A|, 

> 1. 


(3.29) 


From these estimates and (I3.17I) - (|3.19D we get 

I /a,21 < |A| Q_1 , 

and 

i/a, 3 i < (i + i^DiAr- 1 . 

We also conclude from (13.91) and (13.291) that 

I'S'a(-) A)| > |A 2 — Aq| = |A + A 0 ||A — A 0 |. 

Thus from (13.251) we obtain the following estimate, valid for each j = 1, 2, 3: 

1 


(3.30) 


I4 + I <)(! + « I 


c\n 


|A - A 0 ||A + A 0 | 


n— (! +/3)J- 


Now we compute J. We split C\f2 into four regions: 


D 1 := {A € C\fi : |A - A 0 | < |A|, |A - A 0 | < |A + A 0 |}, 

D 2 := {A € C\n : |A - Aol > |A|, |A - A 0 | < |A + A 0 |}, 

D 3 := {A € C\Q : |A + A 0 | > |A|, |A - A 0 | > |A + A 0 |}, 

D a := {A € C\n : |A + A 0 | < |A|, |A - A 0 | > |A + A 0 |}. 
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We split the integral J accordingly 


; = Et 


| A — Ao 11A + Ao 11A | 1 ~ 

D 




3 = 1 


Now we are ready to give uniform (w.r.t. u ) estimates of Jf 


1 


w s t 
1 


1 


1 

< - 


1 


< I 

i "> 


t J |A — A 0 | 3_ “ “ t J |A| 3 -“ ~ t 

Di |A|>1 

1 


< 


1 


^ - t J |A|3-“ “ t / |A|3— a ~ t 

D 2 |A|>2 


, T , 1 f 1 1 

3 “ t J |A|3—« - ~t 

d 3 I. 


1 


\Ji\< 


[ - [ _ 1 

' J | A + A, 

Di 


|A|>2 

1 


|A| 


3—a 


<1 


- < - / --- < - 

0|3— a ~ t J |A|3"« ~ t 

| A |>1 


Note that these estimates are valid only if a < 1. We conclude that 

1 


ifji < a + 


t 


The estimate for I 3 is carried out similarly. First of all, note that from (13.911 


S\\ = d\ 


^|(A 2 - Aq)(A 2 - A 2 )(A 2 - A 2 ) 


E 

.7=0,1,2 


^a + 5a 


(A — A j) (A + A j) 


-4 


5a 


— : 5 aa + • • • + 5j A . 
Further, on C\fi, and using (13.30(1 . 


51a 

<r 1 

Six 

<r 1 

(5a) 2 

~ A — A 0 2 |A + A 0 ’ 

(5a) 2 

~ A — A 0 A + A 0 2 


Six 


(5a) 5 


< 


|A - A 0 ||A + A 0 | 


j = 3,--.,7. 


(3.31) 


(3.32) 


From now the estimates are obtained similarly to the way they were obtained for I 2 ; more 
precisely, by splitting the domain of integration into the subdomains D\,... ,D 4 . We conclude 


ihi<) 


|/(A)||S aa («,A)| < 1 


1 


c\n 


|5 a («,A )| 5 


<1 


t J |A|>1 |A | 3 “ ~ t ' 


(3.33) 


Now we pass to estimating integrals I±, I 2 , and I 3 . The reference domain now is £l\(D e U 
.Bi(O)), described in (13.2411 and Fig. SI First of all, let us denote 


A(!) = A 0 , X {2) = -A 0 , A<3) = Ai, A< 4 > = —Ai, A® = A 2 , A® = -A 2 . 
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Figure 4: The dark region is O\(i4 e URi(0)). The curves composing dD e \B i(0) in I\ (see (13.141) 1 
are the darkest continuous lines. 


We split the integral I\ into the sum of the following integrals 

4 - "" D tS(u,\) 




fc=i 


me* 


dBe^W^BiiO) 


SxW 


-dX 




k=1 


Further, we split Q\(D e U Ri(0)) into four domains (k = 1,... ,4) 

:= {A g n\(Si(0) U D e ) : |A - A^| < |A - A^|, £ ^ k, 1 = 1,. 


••,4}. 


(3.34) 


(3.35) 


Note additionally the essential inclusion, valid for each k = 1,..., 4, 


dB e (X^)\Bi(0) C n (k) . 


(3.36) 


Finally, we put 


and 


k = 1 f h{X)e itS ^ 

2 ’ t J S\(u, A) 

nw 

T -,k , = 1 f f(X)S xx (u,X)e itS ^ 

3 ’ t J {S x (u, A)) 2 

n( fc ) 


(3.37) 


(3.38) 


In order to estimates these integrals and each l \, we follow the approach that we described now. 
For every integral over domain fl^ k \ we use the following polar coordinates at the stationary 
point A( fc ); 

A = A (fc) + pe 1 ^, p € [0, 2ir), 

and where e < p < po{p)- Note that po is always bounded by a fixed constant Sq, that can be 
chosen e.g. equals 4, uniformly on p. 
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Now we start computing. Fix p as previously mentioned. 

1 . On dB p (X^) 0 we have 1 < |A| < |Ao| + 1, so we use the following estimates: 

<R, 


A 4 — 1 

< 1 

A 2 — 1 

A 4 


|A| 


(recall that R = |Ao|, defined in (13.201) '). On the other hand, using (13.91) . 

\Sx\ ~ iym |A - A«||A - A( 2 )||A - A( 3 )||A - A^||A - A( 5 )||A - A( 6 )| > P R , 

l A l 


and from (13.321) . 


so that 


Saa 


1 0 

JrE- 


4 i 


(S x y 5a A-A0-) AJ ’ 
5aa 




< _L + _L 

~ pR? p 2 R 


Using these observations we are able to obtain the following estimate for l\ in (13.341) : 

1 


Wl * 7 


i 


dBR\m)\BR0) 


|5a(A)| 


< 1 R a X e = l^a-1 
~ t eR t 


On the other hand, for (I3.37D - (I3.38I) . 


p < 1 f |/a(A)| 

12 1 ~ t J |5 a («,A)| 

eo -i 

< la + m j ~ j = 1,2,3, 

£ 

and finally, 

x < 1 f |/(A)||5 A a(«,A)| 

3 ~ t J \S x (u, A )| 2 

n(0 

^ 1 J R a pdp 1 J R a pdp 

~ ij pR 2 + 1J p 2 R 

£ £ 

< -i? Q - 1 |lne|. 

~ t 11 

(Recall that eo is a certain fixed independent of u. t, e, and R.) 

2. On dB p ( A^ 2 i) D fl( 2 i, computations are perfectly symmetric with respect to the first case. We 
obtain the following estimates 


I If 


<±R a -\ |/ 2 -’ 2 |<^(l + /3)R“- 2 , 


j = 1,2, 3, and 


I4"’ 2 I < -R a ~ l \ 

3 t 1 


lne|. (3.39) 
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3. Now we deal with the more involved case dB p {X^) n In this region we have |A®| = 1 

and 1 < |A| < 2. Therefore, from the decomposition A = A® + pe we shall use the following 
estimates 

l|A| 2 -l|<P, ||A| 4 -l|<p, |/(A)| < p a+1 . 

Similarly, 


^^|A-A ( 1 ) ||A-A( 2 )||A-A( 3 )||A-A( 4 )||A-A( 5 )||A-A( 6 )|>pl. 2 . 


Finally, since 


we obtain 


S\\ 1 r 




£ 


' — 4 1 


(S'a) 2 5 A L^A-A(i) AJ’ 


5 A a 


(5a ) 2 


< 


1 


1 1 . ^ 
1 4-H — < 


1 


+ 


1 


+ 


1 


pR 2 \ p RJ ~ P R 2 p 2 R 2 pR 3 ‘ 

Using these observations we are able to obtain the following estimates: 

1 e■e a+1 1 s a+1 


|/?| < 


t £R 2 t R 2 ’ 


£0 


iri,3i ^ 1 n | m [ P a P d P (! + 

|J * ls h 1+w J yjp- rie 


3 = 1,2,3, 


£0 


I III < 


0 Q ;+1 


1 1 1 i , . 1 

+ ■VTPPT + —pap < 


IpR 2 p 2 R 2 pR 3 J' 


tR 2 ' 


4. Finally, the case for dB p (X^) n O^ 4 ) is very similar to the previous one. It is not difficult to 
see that the following estimates can be obtained 


i/ 4 | < i/-- 4 | < ^ + 

1 i 1 ~ t R 2 ’ 1 2 1 ~ tR 2 


j 1,2,3, \I 3 ’\< tR2 - 


1 


Adding the four previous estimates, we conclude that 


and for k = 1,..., 4, 


and 


3 =1 




R a + 


R J 


i4’“ ,fc i<Ei 7 ii^ 

3 = 1 


(1 + 


R 2 t 


(i + R a ), 




3 = 1 


(3.40) 


(3.41) 


(3.42) 


Collecting all the previous inequalities, including (13.2311 . (13.311) . (13.331) and (13.401) - (13.421) . we 
get the global estimate (recall that R = |Aq| > 2 and t is w.l.o.g. assumed large) 


2 d a , (1 + 1^1) 


\I\<e Z R a + 

L 

< + ii± 


<r“+l 1 1 

1 + R a ~ l + + —r + R a ~ 2 + - + i?" _1 |lne| 


R 2 R 2 

F a+1 

1 4——^—b R a 1 1 In er | 
R z 


R 
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Figure 5: Setting for Case 2 (in fig., Lemma 13.11 item ED- D e is the dark shadowed region. 


By choosing e := min{^=, ^=}, with 5q > 0 small, we obtain that, if a < 1, the integral I is 
uniformly bounded with respect to R and decreases as 




^ r *- 1 (i + 


tR 


1 + In tR 


(1+ |/3|)| lnt| 


3.4.2 Case 2: u € U 

Recall that in this case all the stationary points lie on the unit circle (see Lemma 1.3,11 items 1, 
2 and 3, and Fig. 

As we mentioned earlier, we take D e to be the union of disks with radius e and with centers in 
the stationary points. We start by estimating R n t in (13.121) . First of all, note that on D £ \B i(0) 
the following estimates are valid: 


Thus, using (13.111) . 


and 


(|A|^-ir< £ Q , |A| ~ 1, |A| 4 -l<e. 


I/(A)| 


^ HA|^-ir(|A| 4 -D ^_ Q+1 


|A| 


cx+4 


\hnt\< J |/(A)|<e Q+3 . 

D £ \B!(0) 


(3.43) 


We pass to the estimate of I ex t introduced in (13.131) . Recall that I ex t is composed of three 
different parts, I \, I 2 and I 3 , see (|3.14|) . (13.151) and (13.161) . 

As in the previous case, we will start by splitting each of the integrals P 2 and I 3 into two 
additional pieces, the first one being an integral over the set of A at a distance greater than some 
fixed value from the unit circle (containing all the stationary points), and the second one over 
the rest of the domain included in C \(D e U B i(0)). More precisely, in this case we set 


C\fi := {A: |A| > 2}, 
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Figure 6 : Case 2 . The white region is C\f2. The light shadowed region (without small disks of 
radius e around the stationary points), is Q\(D e U F?i(0)). 


and we denote as I 2 . Tf, the integrals with same integrands as I 2 and I 3 , respectively, and with 
domain of integration being C\f2 for the plus sign, and £l\(D e U -Bi(O)) for the minus sign (see 
Fig. [ 6 ] and Case 1 for a similar splitting). 

We first treat integrals I 2 + , I 3 ■ For that we use the following estimates valid on C\fi: 


A 4 — 1 

< 1 

A 2 — 1 

A 4 


A 


|A|, 


and from (13.91) . 

\Sx\ > |A| 2 - 

From the decomposition (see previous case) 

c / sr' (S\(u,X) S\(u, A) 

S « ( “’ A)= A- ((wX) + (Ary) 


J=0,l,2 


— 4 


Sx(u,\) 

A 


we also have 


S\\(u,X) 


< — 

~ |A| 3 ' 


(s x (u,x)y 

Finally, using (13.111) and (13.1 71) - (13.191) . and estimates (13.441) . we have that 

I/(A)| < |A|“, |/a,i(A)| + |/a, 2 (A)| + |/a, 3 (A)| < (1 + IflJIAr 1 . 

Thus we obtain the following estimates (recall that a < 1) 


< 


1 


I/a(A)| 


< 


1 


1 c 1 


I4 + I < 


tJc\n\Sx(u,\)\~tJ |A| 3_q ~ t ’ 

c\n 

1 f |/(A)H5aa(u,A)1 ^ 1 f _J_ 

t J \S\(u, A )| 2 ~ t J |A|3-“ ~ t 

C\f2 C\f2 


< 


1 


(3.44) 


(3.45) 

(3.46) 
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Now we pass to estimating integrals I±, I 2 , / 3 . As in the previous case, we denote 


A« = A 0 , A^ = -A 0 , aW = A 1; A< 4 > = -Ai, A< 5 >=A 2 , A< 6 > =-A 2 . 


(Note that they are not necessarily different.) Now we split the integral I\ (see (13.141) 1 into the 
sum of the following six integrals 


h 



E 


k 1 as e (A('=))\B 1 (o) 


f(X)e tS ^ u,x ' ) 

SxW 


dX 


=■ E 4- 


Further, we split £l\D e into six domains, for k = 1,..., 6, 

= {A E fi\(D e U5i(0)) : |A - A (fc) | < |A- A W |, l^k, 1= 1,... ,6}. 


Following the previous case, we define 

k = 1 f f x (X)e tS ^ 

2 t J S\(u, A) 

n( fc ) 

T -,k = 1 f f (X)S\\(u, X)e tS ( u,x ' ) 

3 t J (S\(u, A)) 2 
n( fc ) 

For every integral over domain we perform the change of variables (see previous case for 
details) 

A = A (fc) + pe itf> . 


Now we compute each integral. Fix p > 0. On dB p ( A^')) n we use the following 

estimates: 

|A| 4 -l<p, |A| 2 -l<p, |A|~1, 


and 


On the other hand@ 


l/(A)t<p“ +1 , |/j(A)| < (1 + \B\)p“. 


|Sj| ~ hjIA - A«||A - A«||A - AW||A - a1 4 >||A - A< 5 >||A - A< 6 >| 


Finally, 


5 A a 

(5a) 2 


1 6 1 1 
— \Y _-— -4— 

S\ L A — A(J") A. 
j=i 


7 In the worst case, where three stationary points coincide, this estimate is optimal. Note also that the other 
three stationary points are antipodal, lying at distances ~ 1. 
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Using these observations, and the inclusion dB e (\^)\Bi(0) C we are able to obtain the 
following estimates for each k: 

f l/(A)| 


I if I < 

I 1 I 


as £ (A( fc ))\Bi(0) 


|5a(A)| 


dX 


for j = 1,2,3, 


\I~’ k \ < 

I 2 I 


1 P a+lp 1 

< ££ _ £ _ 

~ t e 3 t ’ 


I/a(A)| 


t J |S a (u,A)| 
n<*0 


< 


£0 


< l(i + IfflK- 1 , 


and 


13 ~ -j- 


1 


< 


|/(A)||5 aa (u,A)| 
\S\(u, A)| 2 

1 f p a+l pdp 1 x 


U 


< -f 

t 


(3.47) 


(3.48) 


(3.49) 


where £q is a certain fixed constant, independent of u and t. 

Gathering estimates (13.43ft . (13.45ft . (j3.46[) and fl3.47ft - fl3.49p . we obtain 

m^e a+3 + ^(l + (l + l/3|)e“ -1 )- 

Finally, for t > 0 large we choose e := min j<5o, ^74 | to obtain that, if a < 1, then 

m<a+ 


1 1 

+ 


< 


(1 + 


t t( a+ 3)/ 4 J ~ t(“+3)/ 4 ' 


3.4.3 Case 3 


For the sake of easiness, we recall the definition of case 3: the region u E C\U and |Ao| = 1 +ui < 2 
(see Lemma EU item 4 and Fig. 0. 

As usual, we take D e to be the union of disks with radius e and with centers in the stationary 
points. For t > 0 large enough, we consider the following two sub-cases: 


0 < co < 


U/4’ 


and u > 


U/4- 


1/ 0 < u < ^74. 

In this subcase we set e = 7^ thus obtaining that oj < e. 

We note that the following estimates are valid on D e \B i(0): 

(IA| 2 - 1)“ < (a; + e) a < e Q , |A| 4 - 1 < a; + e < e. 


Further, it is easy to see that the reasoning used in the previous case applies to this case 
also. Thus we obtain that |7| < —^3-. 
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Figure 7: The setting for Case 3. Note that if e.g. Ao, Ai and A2 are close, and A is close to them, 
then —Ao, —Ai and —A2 are at distance ~ 1 of A. This fact is useful to get a lower estimate on 

5a- 

11/ ca > ^. 

In this case we set e := ^+3 thus obtaining that 

2 

CJ > 4 £-“+3 > 4 £, 


The integral Ii n t from 1)3.1211 is directly estimated as follows: 

\hnt\< [ l/(A)|<e 2 <^y. 

J t~ 

AA-Br(O) 

Similarly, the integrals l/, Z/~ are estimated as in Case 2, u € TJ. 

Now we pass to estimating integrals Ji, 7^”, 7^~. Similarly to the previous case we split 
each of these integrals into a sum of six integrals and perform the change of variables of 
the form A = A^ + pe 1 ^. 

In order to estimate I\ we use the following estimates valid on dB £ (\( k >) n 
(|A| 2 - 1) Q < (w + e) a < u a , |A| 4 - 1 < a; + e < a;, 
l/(A)l < u,°+\ |/a(A)| < (1 + I/3IK, 

ancd 

|5 a | ~ p^|A - A«||A - A^||A - A^||A - A^||A - A^||A - A^| > ew 2 . 

8 The worst case corresponds to uo being close to 1, and |A 2 1 is chosen close to 1 
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This allows us to obtain 


U.i < 

1 11 ~ t £U 2 

1 


< 


tu> 


1—a 

i 


< 

~ 2(1 — a) 

ts a + 3 

1 < 1 


t 1 4 t 4 


1 —a r-u ai-|-3 


In order to estimate I 2 ’ , / 3 : we use the following estimates on 8B p { A^) n 


(|A| 2 - l) c 


< 


u a , if 4p < u) 


, |A| 4 - 1 


p a , if 4/9 > UJ 

For the term S\, we use the fact that in the considered region, 


ui, if 4 p < uj 
p, if 4 p > uj 


s x ~ -L|A - AW||A - A^||A - A®||A - A^||A - A^||A - A^|, 
|A| 


to conclude that 


\Sx\ 


> 

r^i 


puj 2 , if 4p < u 
p 3 , if 4 /j > uj 


Finally, following the same approach as before, we have 

5aa 


(5a) 2 

In this manner we get the following 


< f pfa’ if A P < w 


jt, if 4p > u 


-,fc, / 1 [ P a+1 dp , 1 f u a pdp 


l^’l < 


i r p a+l dp i r 

t J p 3 + t J 

u/A £ 


< 


< 


pu: 2 ~ tuj l ~ a ~ ^ ' 


and 


£0 


uj/A 


-,k\ [ P a+2 dp | 1 f u a+1 pdp < | lne| ^ [ In t| 


I /“’"I < 

I O I 


* + t 


< 


uj /4 


9 2 w 2 ~ to; 1 - 0 ~ ' 


Thus we have proved that in this case 


\I\ 


< IM 

rsj q+3 "> 

t 4 


as desired. 
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3.5 Proof of the estimate for small t 

3 

We will suppose that t < e 1 ~ a (which implies, in particular, that t < 1). Note that in this case 

-j- > e and thus for all uj > — r we have that is a decreasing function of u; and thus the 

£3 £3 u 

following inequality is true: ^ for all uj satisfying the above-mentioned condition]^ 

U 3 i^T- 

We consider the following cases for the values of parameter u. 

1. All the stationary points lie in ^(O). 

2 . Stationary points Aq, — Aq lie outside the ball ^(O). 


3.5.1 Case 1 

In this case we set D r 
two parts: 


Evidently, we have 


B r ( 0) with r = -573 (note that r > 4) and we split the integral I into 


/ + /’ 


I — / T / •— lint T lext- 

D r \Bi(0) C \D r 


Iint\ < r a+2 . 


Integrating I ext by parts, we represent it as a sum of three integrals (see page 12): 


lext = 4+4 + 4 - 


On C \D r the following is true: 

\Sx\ ~ |A| 2 . 

Indeed, 


A 


> 


1 




|A — Xj\ > |A| - |Aj| > |A| -2 > |A| 


|A[ 

2 


]A| 

2 ' 


In a similar way we have that 

S\ a ^ 1 
(Sa) 2 ~ |A|3 


This allows us to obtain the following estimates: 


141 


< — 

1 5 

tr 


\h\ 

141 



1 

tr 1_ "’ 


1 

tr l ~ a 


Finally, taking into account that r ~ 773 , we get that |/| < . 

t ' ,— n — 


9 This assumption is used to treat In to arising in (I 3 . 50 I) - (I 3 . 51 II . It can be replaced by the assumption that t < 1 . 
Then, we can note, for example, that for all £ > 0 and to > 1 we have that lnui < to 6 . Then at the end we obtain 
the following estimate: |/| < Q+ 1 e +'i ■ Then, if e < 1 ^, we also have that |/| < -777-. 
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3.5.2 Case 2 


In this case we consider two subcases 

1. 1 < u < -jYj 3 In this case we take r = and set D r = B r { 0). Note that on C \D r the 
following estimates hold: 


|a ~ > i _ M > i _ 1 + uj > I 

|A| - |A| - r - 2’ 

|A — Xj\ > |A| - |Aj| > |A| - (1 + u) > |A| - - > y-. 


Thus, in this case the estimate can be carried out by using the reasoning of the previous 
case. 



In this case we choose r = min 


l 

mr 




and we set D r to be the union of B 1+ 1 (0) and 


.B r (±Ao). Note that if -^+4—— > \ (and hence r = 4), we have that u a < 

t c jl 

we split I into Ii n t and I ex t as above. First of all, note that 


l 

a + 2 ■ 


Further, 


L int | 


< i + w 


a 2 


< 


1 


a+2 • 

t 3 


Now we introduce the following notation: B p = Bi +p (0) U -B p (±Ao). Note that 

C \D r = (C \B 1/4 ) U (S 1/4 (±A 0 )\5r(±A 0 )). 

Each of the integrals /*,, k = 1,2, 3, is split into the sum of two integrals I^ + I p . where 
I^ is the integral over the set S 1 / 4 (±Ao)\S r (iAo) and I£ is the integral over C\Bi/ 4 . 0 

Note that on dB p with p > | and on dB p (± Ao) with r < p < the following estimates 
are true: 

IA ~ A s t| . . _ I Agf | ^ 

|A| - |A| ~ ’ 

where A s t = ±Ai or A s t = ±A 2 - Thus we have that on dB p with p > \ and on dB p (± Ao) 
with r < p < A the following holds: |S,\| > |A — Ao11A + Ao | • Using this fact we can obtain 
the following estimates valid on dB p with p>\ and on dB p (± Ao) with r < p < 


|SaI> 


up, 4 p < u, 
p 2 , 4 p > u, 



4 p < u, 

U)p z ’ r ’ 

js-. 4/5 >U. 


Using these estimates, we are brought to the following estimates on the parts of I ex f- 


iyt< 


u 


tu 


1 ^ 1 

tu 1 - 01 - t +r ’ 


u a r 1^1 
lur ~ tu l ~ a ~ ’ 


10 To be more precise, if r < then + is the integral over dD r n Bi/ 4 (±Ao)\Br(±Ao), i.e. over <9B r (± Ao) and 
I j" is the integral over dD r i.e. over ctB 1+1 / 4 (0). Note that if r > 4, then /*, = l£ and I p are formally 

not defined in this case, k = 1,2, 3. 
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Figure 8 : The setting for Case 2, subcase 2. The set B r is the shadowed region. The black 
circles form dB p (± Ao) with r < p < the dashed circles are dB^/ 4 (±Ao) and the dotted circles 
form dB p with p>\. 


<*>/4 oo 

f u a ~ l pdp r p a ~ 1 pdp < 1 _j_ 

] tup + J t{? ~tu'-«~ t *¥' 

1/4 u/A 





r 


1 

— a+2 ) 

t 3 


u/A 


Wl S 


u a pdp j' p a pdp < lncu ^ | In 11 


1/4 

1/4 


tup 1 


u/A 


tp A 


< 

tU 1 “ Q ~ ’ 


\Is\< 


u a pdp | In r | ^ | lnt | 


< 


tup 2 ~ tu l ~ a ~ ^ ' 


(3.50) 


(3.51) 


4 Bilinear Estimates 

The purpose of this section is to show bilinear estimates for the NV equation for negative energy, 
always taking into account the size of the fixed energy. 

In this section we use the following notations. We denote by (/) the japanese bracket: 

(/) ; =(i + l/l 2 ) 1/2 ; 

J-[u] and u both denote the Fourier transform of u in (t, x, y)\ we also have AAB:= ruin (A. B) 
and Ay B := max (4, B). Variables N , N, N, L, L, L of this section are dyadic , i.e. their range 
is {2 k ,k G N}. 
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Now we introduce the associated X s,b spaces [5] for the NV dynamics. Let (p E Co°(M; [0,1]) 
be a cutoff function such that 


ip(s) = 0 for |s| > 1 , 


<p(s) = 1 for 



We define 

‘ P( s ) ■= v( s ) ~ v( 2s )- 

We introduce the frequency projection operators at a dyadic frequency N > 1, as follows: 

<Pn(s) ■= ipiN^s), 

and for N = 1 , 

¥>i(«) : = £(«)• 

Using these multipliers, we have for the Fourier transform in the £ variable, 

p NU ■.=. 


Recall the definition of the phase S(u,f) in (12.21) . We define 

:= (£ 3 + f) (i-p) • 

Note that the above expression is real-valued. In order to perform some Fourier analysis, we 
need w in terms of real-valued coordinates. Put £ = £i + i£ 2 . We have 

M£,£) := w(t,t;E) = 2(£ 3 - 3£r£ 2 ) (l + . (4.1) 

Compared with the work of Angelopoulus T| , the NV± symbol has a nontrivial component 
which becomes important at low frequencies 11 l. We also define 

o-(u£ 1,£2) := t - ui(£i,£ 2 ). 

We introduce then 

Q l u := E-\ip L {\E\-^\)n^]{r^)}. 

Finally, for a fixed energy E, we say that u = u(t, x, y) E X^ b for s, b E M if u E L 2 (M 3 ) and its 
Fourier transform u satisfies the integral condition 


(o-) 26 (|£|) 2 s |u(r,£i,£ 2 )| 2 dTd£id£2 < +00. 


(4.2) 


Proposition 4.1. Assume that E < 0. Then we have for e > 0 small and s > 
\\d z (vw)\\ x s,-i/2-2e < |.E | (3 ~ 4s)/8 \\ v \\ x ^ 1/2+s \\w\\ x ^ 1/2+e , 
for all v , w such that the right hand side makes sense. 


(4.3) 


1 Although compared with the usual KP symbols, the NV symbols are still bounded at the origin. 
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Proof. We follow closely the ideas from [23j, using a modified version of the original ideas by 
Bourgain [5] and Kenig, Ponce and Vega m- As usual, by duality we are lead to prove that 

J := K[t,t,£,£\u(t, 6,6)#,f:I,6)d£i 2 d £ 2 dr df 

Jr 6 

satisfies 

l+l 6 II^IIl 2 II^IIl 2 II^IIl 2 ; 

for any u,v,w € L 2 (M 3 ), and where for j = 1,2 , a := cr(f,£), 

f := f — r, i ■= i — £, a := a — a. (4-4) 

Here the kernel K = K[t, f, £, £] is explicitly given by 

k : = i^i+i2ii^i+6ii - &+6i- 1 (^}- 1/2+2 6^)- 1/2 - E (^>- 1/2 - £ (iei} s (ier s (iir s . 

A further simplification leads to 

I< = |f|W- 1 '' 2+2 '(<r)- 1 /2-=( # )-i/2-«(||| ) . ( | {r . ( |j|)-». (4.5) 


Now we use dyadic decompositions to split J into several pieces. We put 

J = ^N,N,N' 

N,N,N 


where 

J nnn := / 4l[r,f,C,|]PAru(r,6,6)+7v v ( f )6>6)+ , i v' u; ('r>li>6)d6^lirf6rf6dr(if. (4.6) 

li 6 

Estimate for low-low to low frequencies. This is the simplest case. Here we have N ~ N ~ N, 
where N ~ 1. From (14.61) we have to estimate the quantity 


E 


J 


N,N,N‘ 


N~N~N~1 


Note that we also have 


|£| ~ ||| ~ ||| ~ | E \ 1//2 . Therefore, using 63D we have 

\K\ < | E\^- s ^ 2 (a)- 1/2+2£ (a)- 1/2 - £ (a)- 1/2 ~ e 
< |^;|( 1 - s )/ 2 (cr)- 1 / 2 _e ((j) _1 / 2-e . 


Therefore, from (|4.6I) . Plancherel, and using Cauchy-Schwarz, 

\J N ^\ 6 l^l (1 - s)/2 E-^iar^P^u] | l4 7- l [{a)- 1/2 - £ P^v] 


L 4 


I PatW 


n w Wl 2 ■ 


From (12.71 ) and the definition of X s,b space norm, we obtain for e > 0 small, 


P~ l [{a)- l / 2 - e P^u] 


L 4 


Z 4 


C/(t)t/(-t)^- 1 [(a)- 1 / 2 - £ P^i] 

< \E\~ 1 / 1&+ ||t/(t)C/(-t)^- 1 [( C r)- 1 / 2 - £ i^u] 

< \E\~ 1 / ie+ \\{a)^ + {a)~ l ^ 2 ~ £ Pnu\\ L 2 

< \E\- l l 1&+ \\P N u\\ L 2 . 
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We conclude that 


E 


^N,N,N ~ 


\E\ 


(3“ —4s)/8 1 


“IIl 2 IfIIl 2 IFIIl 2 - 


N^N^N^l 


Estimate for low-high to high frequencies. The worst case here corresponds to N <C N and 
N ~ N, where iV 3> 1. From (14.6p we have to estimate the quantity 


E 


J 


N,N,N' 


jv<siv, !V~jV»l 


Unfortunately, a crude estimate for K in (14.51) shows that it is not possible to counter balance 
the term |C| ~ \E\ l ^ 2 N. For this reason a dyadic decomposition on the modulation variables 
is necessary. Using the frequency localization operator Ql, Q^ and Q^ on dyadic shells a ~ 
\E\ 3 / 2 L, a ~ \E\ 3 / 2 L and so on we have 


J; 


N,N,N 


N,N,N’ 

L,L,L 


where 

:= J K [T>f,Z,£]F[PNQLu]{T,Z U S2) X 

xE[P^Q l v ](f, li ,£ 2 )J : [PnQi w \(t, £1, £2) df 1 df 1 df'idf'zdr df. ( 4 . 7 ) 

We readily have in the considered region 

\K[t,t,£,£\\ < \E\-( s+7)/a NN- s L- 1/2+2£ L- 1/2 - £ L~ 1/2 - £ , 
and using Cauchy-Schwarz, the fact that 


/ P[P n Qlu](t, Ci ,&)F{Pf l QLv]{f, Ci, £2) 

J t,£,t f 

represents a convolution in Fourier variables, and Plancherel, we get 

\ J n L P L n\ ~ |^|- {s+7)/4 iVAr- s Z- 1/2 + 2£ L- 1/2 - £ L- 1/2 - E ||P J vQ L ir P^Q l v\\ L 2\\P^Q- l w\\ L 2. (4.8) 

Now we use the following estimate (see [23] for a similar statement), valid under the assumptions 
that we work with low-high to high frequencies: 

II PnQlu P„Q l w\\ L 2 < \E\ 5 / 4 N 1 / 2 N~ 1 L 1 / 2 L 1 / 2 \\PnQlu\\l2 \\P^Qiw\\ L 2. (4.9) 

This estimate is proved several lines below. For now we assume the validity of this estimate and 
we continue with the estimation of (14.81) . Note that (14.91) allows to cancel out the bad frequency 
N in (14.81) . We get 

\ J n L P L n\ < | J E7|- (2 + s)/4 A^ 1/2 - s Z- 1/2 + 2 " J L--Z--|| J P iV ^^||^|| J P Z vQz / ^ll^ll^vQz^ll^- 

Adding up on N (s > |), L, L and L we obtain 

^N,N,N ~ l-^l ^ + ^ IMIl 2 y! II^JV^IIl 2 II-Pjv'HI-L 2 
1V<JV, 1V~JV> 1 N~N 

< \E\^ 2+s)/i \\u\\ L 2 \\v\\ L 2 \\w\\ L 2 . 

12 The power 3/2 of the energy in front of the modulation variable a is dictated by the time scaling. 
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Let us prove (14.91) . Following [24) . and using the Plancherel’s identity, together with Young’s 
inequality for convolutions, we have 


\\PnQlu PjyQiw \\ L 2 ~ \\E[PnQlu\ * F[Pf]Qiw \\\ L 2 

< sup (measA E (f,£)) 1 / 2 \\P N QLu\\ L 2\\PfiQ L w\\ L 2, (4.10) 


where A E (f,£) is the set (see (14.41) 1 

Ae{ f,i) ■■= {(r,0 € M 3 : M ~ \E\ 3 / 2 L, |<j| ~ \E\ 3 ' 2 L, |£| ~ {E^N, ||| ~ \E\ l / 2 N). 
The measure of this set can be estimated as follows: 


meas A e (t , £) < \E\ 3 ^ 2 (L A L) meas B e (t, £), 


(4.11) 


where 

Be{ f,i) := {(r,0 €R 3 : |d + tf[66l < |£| 3/2 (LV L), |£| ~ \E\^ 2 N, ||| ~ {E^N}, 
and E[ is the standard resonance function (see (14.11) 1 


H[£, |] := w(£, £) - w{£, £) - w{£ - 6 £ - 0- 

We use now the idea from [241 Lemma 3.8]: we estimate the measure of B E by finding lower 
bounds on the derivatives of the function £ H > #[£,£] in the considered region. 

As in pp, we have to distinguish between two cases: |6 — 6| 16-61 rs j \E\ l/2 N ~ I E\ l / 2 N 

and |6 — 61 ^ |6 — 61 (the remaining case is identical). For the first case we have 


d i2 H[ 61 ] = - 2 % 


(£? - 366 2 ) ( 


1 + 


3|6 


= -12 


v 6 + 6 

+((6 - 6) 3 - 3(6 - 6)(6 - 6) 2 ) ( i + 
3|6 


3|6 


- 66 i + 


6 2 + 6 


+ (6 - 6)(6 — 6 ) i + 


(6 - 6) 2 + (6 - 6) 2 
316 


(6 - 6) 2 + (6 - 6) 2 


1666(6-36) , 16(6-6)(6-6)((6-6) 2 -3(6-6) 2 ) 


( 6 2 + 6) 2 


+ 


((6 - 6) 2 + (6 - 6) 2 ) 2 


Using the estimate 


ab 


a 2 + b 2 


~ 2 ’ 


and similar other estimates for the fractional terms appearing from the fact that we work with 
nonzero energies, and valid for all (a, b) E R 2 (the limit at the origin is not well-defined, but the 
functions are always bounded), we get 


\dt a H[£,i]\~\E\N 2 
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In the second case, we have no problems since 




= -2d, 


& 


(£1 - 3 6 <£l) 1 + 


3|^| 


= -6 


« + « 

+((li - 6) 3 - 3(6 - 6X6 - 6) 2 ) (i + 

-((6 -fi) 2 -(6-6) 2 )(l + 


3\E\ 


(6 - 6) 2 + (6 - 6 ) 


3|£| 


( 6 - 6) 2 + ( 6 - 6) 2 


meM - 3 e 2 2 ) + mdi - 6) 2 «6 - 6) 2 - 3(6 - 6) 2 ) 


( 6 2 + 6 2 ) 2 


((6 - 6) 2 + (6 - 6) 2 ) 2 



so that 

|%tf[ 6 £]M£|iv 2 . 

We conclude (see (23J for example) that 


measB E (f,i) < \E\NN~ 2 (L V L). 


Finally, from (14.1111 and (14.1011 we conclude. 

Estimate for high-high to low Jhh^l, and for high-high to high frequencies Jhh^h- This is 
the difficult part of the proof, because for obtaining m with s > 4 we do not have the 
corresponding Carbery-Kenig-Ziesler [7] result. Instead, we will prove that the corresponding 
smoothing estimate Lemma 12.11 suffices for the case of negative energies. 

We prove the most difficult case, the one for high-high to high frequencies (see below for a 
comment on the case high-high to low). Here we have N ~ N ~ N, where N 1. From m 
we have to estimate the quantity 

Jn,n,n- (4-12) 

AT~iV~iV»l 

Note that we also have |£| ~ |£| ~ |£| ~ \E\ l / 2 N. Therefore, using (14.51) we have 
\K\ < \E\~ S / 2 \Z\ lV- s (d)- 1/2+2e ^) _1/2_£ (d) _1/2_£ 


Therefore, from (14.611 

\ j n,n,n\ £ | E |- s / 2 +V 4+ jv 1 / 2+ “ s |3^| 1/4 ~J r-1 [(cj)- 1/2 - £ P/vu] 


L 4 


|9 2 | 1/4 ~ J r_ 1 [(d') _ 1 / 2 _e P^u] 


L 4 


WIl 2 


From (12.811 we get 


I J 


N,N,N\ 


< \E\~ s / 2+1 / a+ N 1/2+ ~ s F- l [{a)- 1/2 - £ P N u\ 


i- x 


.Y u ’2 


X 


F-'m-^Pxv} 


x0.h~ " W " L2 

< |F1|-' s/2+1/4+ “ 0+ AC 1/2+_s ||Pjvw||l 2 IIT’tv^IIl 2 IIHIl 2 
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Adding on N ~ N, we conclude. 


Finally, some words about the case high-high to low frequencies. In this regime one has 
N ~ N N. Note that we also have |£| ~ |£| ~ |I?| 1 / 2 ./V, and |£| ~ \E\ l ! 2 N. Now it is enough 
to consider the following estimate: 

mi<K 

N S N S ~ N s ’ 

Therefore, using (14.51) we have 


K | < |£?r a/ 2 |^|iV-*(cr)- 1 / 2 + 2 e <<T)- 1 / 2 - 6 (CT )- 1 / 2 - 6 

< | S |-*/2+l/4+ ]V l/2+- S | ? |l/4- {(j) -1/2- £ |||1/4- ^-1/2-^ 


and the rest of the proof is similar to the previous case. 


□ 


5 Local well-posedness 


In this section we prove Theorem ll.il Using the standard Bourgain’s method, we will show the 
following 

Theorem 5.1. Fix E < 0. The Cauchy problem for (II.II) - (II. 21 ) is locally well-posed in LP(M 2 ) 
for s > Moreover, the existence lifetime Te > 0 of a solution v(t) with initial data vq satisfies 

Te\\vo\\%s > IF^- 3 )/ 8 , (5.1) 

for some positive exponent a. Finally, one has for all i) £ [0,1] standard cut-off supported in the 
interval [—2, 2], rj = 1 on [—1,1], and s > \, and T < Te, 

\\v{t/T)v(t)\\ i +£ < ||i7o||jf- (5.2) 

X E 

Remark 5.1. Note that in this result we give explicit dependence on the energy E, for further 
developments. 


Proof. The proof is standard, we only check the main lines of the proof. Assume that vq £ Ff s . 
From Annua) we have the local Duhamel’s formula for f £ [0,1], 

rj T (t)v(t) = rfr(t)T[v](t) := r] T (t)U(t)v 0 + 2 rj T {t) f U(t- t')[d z (vw) + d~(vw)]dt ', 

Jo 

where i] = r){t) € [0,1] is a smooth bump function with rj(t) = 1 for t £ [—1,1], and rj(t) = 0 for 
\t\ > 2, and rjr(t ) := rj(t/T). From this identity, the standard linear estimates (see e.g. Lemma 
3.1, Lemma 3.2 of [9]) and (14.31) we have, for any s > i, 


\\r]Tv\\ x s,i/ 2 +e < \\vo\\h° + C\\r) T d z (vw)\\ x s,-i/ 2 +e 

< ||uo||^ + CT £ \\rj T d z (vw)\\ x s,-i/2+2e 

T £ 

< llwollil- + C |£|(4 s-3)/8 hTvf x s,i,2+e 


Now we fix any time T ~ 


( \E\ ( 4s — 3)/8 

V 8C||«o|| h . 


1/e 

, and the ball 


B:={ve X s f} /2+£ : \\rj T v\\ x s,i /2 +e < 2||u 0 ||h s }- 
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For v € B, one has 


ll^rTHII^i/a+e < 2 ||uo||i? s - 

The contraction property is proved in a similar fashion. The proof is complete. □ 

A Proof of (11.41) 

We have 

(d? + &l)v = ^(9* - idyfv + ^(0* + id y fv 

= id y ) 2 (d x v - idyV ) + ^( 0 * + id y ) 2 (d x v + id y v) 

O O 

= \{dx - id y )(d 2 v - 2id 2 xy v - d 2 v) + ^( 0 * + id y )(d 2 v + 2 id 2 y v - d 2 v) 

= \{d's V ~ ‘ lld L y V ~ dlyyV - i{dl xy V - 2 idlyyV ~ ^V)) 

+^(d x v + 2i 9 xxy v - dlyyV + i{dl xy v + 2 idl yy v - d$v)) 

= \dx{dlv - Sd 2 v). 

Similarly, if w = w± + iw 2 = (w±,w 2 ), 

d z w + d z w = i( d x w - id y w ) + ^(d x w + id y w) 

= d x wi + d y w 2 
= V • w. 

Similarly, 

2 d z (vw) + 2 d z (vw) = 4 Re 0 Z (rm;) 

= 2 Re( 02 ; — id y )(v(w 1 + iw 2 )) 

= 2d x (vwi) + 2 d y (vw 2 ) 

= 2X7.(vw). 

On the other hand, m can be written as 

d y w\ + d x w 2 = 3 d y v, d y w 2 — d x w\ = 3d x v, 

Note that one can recover w\ and w 2 from v via the identities 

Auq = 3 (d 2 v — d 2 v ), AtC 2 = 6d xy v. 

Therefore equation (ED reads now 

d t v = 2[d x (d x v — 3 d 2 v) + X7.(vw ) — EX7.w\, 

as desired. 
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B High-energy limit of NV equations 

Consider the Novikov-Veselov system of equations written in (x, y ) variables: 


d t v = 2d xxx v - 6 d xyy v + 2{d x {ywi) + d y {yw 2 )\ - 2E(d x w\ + d y w 2 ), 

< d x wi - d y w 2 = -3d x v, (B.l) 

d x w 2 + d y w\ = 3 d y v. 

(see also <HZD)- 

Put E = ±k 2 and dilate the second variable: y = kY. Then system (IB. ID becomes 


dfV = 2d xxx v — 6 —zd X YYV + 2d x {yw \) + 2—dy(vw 2 ) — 2Ed x w\ — 2 E—dyw 2 , (B.2a) 

K Z K K 

d x w\ — —dyw 2 = —3d x v, (B.2b) 

K 

d x W 2 + —dyw i = 3 —dyv. (B.2c) 

K K 

Formally, let us search the solution of (IB. 2D in the following form: 


v = v 0 (x,Y), (B.3) 

wi = =f3k 2 +w®(x, Y) +w\{x,Y) — +wf(x,Y)-^, (B.4) 

K Av 

w 2 = w°(x,Y) +w\{x,Y)~ + u|(x,y)4r- (B.5) 

K K- 

We will also suppose that 

Wj(x,Y) —*■ 0 as x —A oo, j = 1,2, A: = 0,1,2. (B. 6 ) 


Formulas (IB.4p . (IB.5D are approximations, up to order —2, of asymptotic expansions with respect 
to large k of solutions to equations (IB.2bp ~ dB.2cP such that w\ ± 3k 2 , w 2 are localized functions. 

Equating terms of order 0 in n in equations dB.2bp . dB.2cD and taking into account condition 
(EH), we obtain 

w® = —3no, W 2 = 0. (B.7) 

Equating terms of order — 1 in k in equations dB.2bD , (|B.2cD and taking into account condition 
EH, we obtain 

=0, w\ = Gd^dy'VQ. (B. 8 ) 


Finally, equating terms of order —2 in k in equations ()B.2bp . ()B.2cD and taking into account 
condition dB. 6 D . we get 


Thus we have 


wl = 6(d x 1 ) 2 dyV 0 , w 2 = °- 

(B.9) 

v = v 0 (x,Y), 

(B.lOa) 

w\ = =f3k 2 — 3v + &(d~ l ) 2 dyv o4r, 

(B.lOb) 

w 2 = Qd x l dyv o-. 

(B.lOc) 


Now let us insert this ansatz into equation dB.2aD . It is easy to verify that the coefficient of the 
term of order 2 in k is equal to zero, and there are no terms of order 1. Finally, equating terms 
of order 0 in k, we obtain the following equation for uq: 


d t v 0 = 2d xxx v 0 - I2vd x v 0 =f 24d x 1 dyV 0 . 


(B.ll) 
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(B.12) 


Note that if vq( x,Y,t) is a solution of equation (IB. Ill) , then 

u(x,y,t ) = -v 0 {-x,2y, ^t) 

is a solution of KP in its classical form 

dtu + 6 ud x u + d xxx u = ±3d~ 1 dyU. (B.13) 

Thus, at high energy E limit, after an appropriate dilation of y variable, equation NV+ 
becomes KPI and equation NV— becomes KPII. 


References 

[1] Angelopoulos Y. Well-posedness and ill-posedness results for the Novikov- Veselov equation, 
preprint ArXiv:1307.4110. 

[2] de Bouard A., Martel Y. Nonexistence of L 2 -compact solutions of the Kadomtsev- 
Petviashvili II equation. Math. Ann. 328, 525-544 (2004). 

[3] de Bouard A., Saut J.-C. Solitary waves of generalized Kadomtsev-Petviashvili equations. 
Ann. Inst. Henri Poincare, Analyse Non Lineaire. 14(2), 211-236 (1997). 

[4] de Bouard A., Saut J.-C. Symmetries and decay of the generalized Kadomtsev-Petviashvili 
solitary waves. SIAM J. Math. Anal. 28(5), 1064-1085 (1997). 

[5] Bourgain J. On the Cauchy problem for the Kadomtsev-Petviashvili equations, Georn. Funct. 
Anal, vol 3 no. 4 (1993) pp. 315-341. 

[6] Calderon, A. P. On an inverse boundary value problem. Seminar on Numerical Analysis and 
its Applications to Continuum Physics. Soc. Brasileira de Matematica (1980), pp. 65-73. 

[7] Carbery A., Kenig C.E., and Ziesler, S. Restriction for homogeneous polynomial surfaces 
in R 3 . Trans. Amer. Math. Soc. 365 (2013), no. 5, 2367-2407. 

[8] Croke. R.; Mueller, J.L.; Music, M.; Perry, P.; Siltanen, S.; and Stahel, A., The Novikov- 
Veselov Equation: Theory and Computation, preprint arXiv: 1312.5427. 

[9] Ginibre, J. Le probleme de Cauchy pour des EDP semi-lineaires periodiques en variables 
d’espace, Seminaire N. Bourbaki, 1994-1995, exp. no 796, pp. 163-187. 

[10] Grinevich P.G. Rational solitons of the Veselov-Novikov equation are reflectionless poten¬ 
tials at fixed energy, Teoret. Mat. Fiz. 69(2), 307-310 (1986), translation in Theor. Math. 
Phys. 69, 1170-1172 (1986). 

[Gr2] Grinevich P.G. Scattering transformation at fixed non-zero energy for the two-dimensional 
Schrodinger operator with potential decaying at infinity. Russ. Math. Surv. 55(6), 1015- 
1083 (2000) 

[11] Kazeykina A.V., Novikov R.G. A large time asymptotics for transparent potentials for 
the Novikov-Veselov equation at positive energy. J. Nonlinear Math. Phys. 18(3), 377-400 
( 2011 ). 

[12] Kazeykina A.V., Novikov R.G. Large time asymptotics for the Grinevich-Zakharov poten¬ 
tials. Bulletin des Sciences Mathematiques. 135, 374-382 (2011) 


38 




[13] Kazeykina A.V.: A large time asymptotics for the solution of the Cauchy problem for 
the Novikov- Veselov equation at negative energy with non-singular scattering data. Inverse 
Problems, 28(5), 055017 (2012). 

[14] Kazeykina A.V. Absence of solitons with sufficient algebraic localization for the Novikov- 
Veselov equation at nonzero energy. Funct. Anal. Appl., 48(1), 24-35 (2014) 

[15] Kenig C.E, Ponce G., and Vega L., Well-posedness of the initial value problem for the 
Korteweg-deVries equation , J. Amer. Math. Soc., 4 (1991), pp. 323-347. 

[16] Kenig, C.E., Ponce, G., and Vega, L. Oscillatory Integrals and Regularity of Dispersive 
Equations, Indiana Univ. Math. J. 40 No. 1 (1991), 33-69. 

[17] Kenig, C.E., Ponce G., and Vega L., A bilinear estimate with applications to the KdV 
equation, J. Amer. Math. Soc. Volume 9, Number 2, April 1996, pp. 573-603. 

[18] Kiselev O.M. Asymptotics of solutions of higher-dimensional integrable equations and their 
perturbations J. Math. Sci. 138(6), 6067-6230 (2006). 

[19] Konopelchenko, B. G. Introduction to multidimensional integrable equations. The Inverse 
Scattering Transform in 2 + 1 dimensions. Plenum Press, New York, 1992. x+292 pp. ISBN: 
0-306-44220-5. 

[20] Lax, P.D. Integrals of nonlinear equations of evolution and solitary waves. Comm. Pure 
Appl. Math. 21 (1968), 467-490. 

[21] Linares, F., Pastor, A., Well-posedness for the two-dimensional modified Zakharov- 
Kuznetsov equation, SIAM J. Math. Anal. Vol 41, No. 4, pp 1323-1339. 

[22] Manakov, S.V. The inverse scattering method and two-dimensional evolution equations. 
Uspekhi Mat. Nauk 31(5) (1976), 245-246. (in Russian). 

[23] Manakov S. V., Santini P. M., Takchtadzhyan L. A. An asymptotic behavior of the solutions 
of the Kadomtsev-Petviashvili equations Phys. Lett. A. 75, 451-454 (1980). 

[24] Molinet L., Pilod D. Bilinear Strichartz estimates for the Zakharov-Kuznetsov equation and 
applications, preprint arXiv:1302.2933, to appear in Ann. IHP Analyse Nonlineaire. 

[25] Murata, M. Structure of positive solutions to (—A + V)tp = 0 in R n . Duke Math. J. 53 
(1986), no. 4, 869-943. 

[26] Nachman, Adrian I. Global uniqueness for a two-dimensional inverse boundary value prob¬ 
lem. Ann. of Math. (2) 143 (1996), no. 1, 71-96. 

[27] Novikov, R. G. The inverse scattering problem on a fixed energy level for the two- 
dimensional Schrodinger operator. J. Funct. Anal. 103 (1992), no. 2, 409-463. 

[28] Novikov, S.P., Veselov, A.P. (1984), Finite-zone, two-dimensional, potential Schrodinger 
operators. Explicit formida and evolutions equations, Sov. Math. Dokl. 30: 588-591. 

[29] Perry, P. Miura maps and inverse scattering for the Novikov-Veselov equation. Preprint, 
ar Xiv.math/1201.2385 [math. AP], 

[30] Saut, J.-C., Remarks on the generalized Kadomtsev-Petviashvili equations, Indiana Univ. 
Math. Journal, vol. 42, No. 3 (1993), pp. 1011-1026. 


39 



[31] Stein, E.M.; and Weiss, G., Introduction to Fourier analysis on Euclidean spaces, Princeton 
University Press (November 1, 1971). 

[32] Sylvester J., and Uhlmann G., A Global Uniqueness Theorem for an Inverse Boundary 
Value Problem , Annals of Mathematics Second Series, Vol. 125, No. 1 (Jan., 1987), pp. 
153-169. 

[33] Taimanov, I. A., and Tsarev, S. P., Blowing up solutions of the Veselov-Novikov equation, 
(Russian) Dokl. Akad. Nauk 420 (2008), no. 6, 744-745; translation in Dokl. Math. 77 
(2008), no. 3, 467-468. 


40 



